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We study the kinetics of the first order phase separation transition in boson-fermion cold-atom 
mixtures. At sufiiciently low temperatures such a transition is driven by quantum fiuctuations 
responsible for the formation of critical nuclei of a stable phase. Based on a microscopic description 
of interacting boson-fermion mixtures we derive an effective action for the critical droplet and 
obtain an asymptotic expression for the nucleation rate in the vicinity of the phase transition and 
near the spinodal instability of the mixed phase. We also discuss effects of dissipation which play a 
dominant role close to the transition point, and identify the regimes where quantum nucleation can 
be experimentally observed in cold-atom systems. 
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Kinetics of the first order phase transitions at ultralow 
temperatures has received considerable attention in con- 
nection with several problems ranging from the decay of 
false vacuum in field theoretical models of the early Uni- 
verse [l| to the phase separation in ■^He-^He mixtures 
In the latter case it has been argued that below a certain 
temperature (of the order of a few tens of mK) formation 
of pure '^He phase from a supersaturated mixture occurs 
via the process of quantum nucleation, where critical nu- 
clei overcome the potential barrier (due to the surface 
energy between the two phases) by means of quantum 
tunneling. It has been predicted theoretically that 
the rate of such nucleation behaves as exp (— C/A/i''/^), 
where A/i is the degree of supersaturation, i.e., the dif- 
ference in chemical potentials of the two phases, and 
C is related to the coefficients of the phenomenological 
Ginzburg-Landau expansion of the free energy near the 
point of the phase transition. Measurements carried out 
by several groups seem to confirm that at sufficiently 
low temperatures the kinetics of the phase separation 
in ■'He-^He becomes temperature independent; however, 
they have been unable to verify the expected dependence 
of the nucleation rate on the systems' parameters (i.e., 
A/i, etc) - partly due to the poor knowledge of micro- 
scopic interactions between particles in such a strongly 
correlated system. 

We argue that contemporary cold atom systems pro- 
vide a perfect setup for studying and observing the ki- 
netics of such a phase separation transition. Atomic 
mixtures, such as boson-fermion mixture, are com- 
monly realized in experiments on sympathetic cooling, 
where one of the species (typically bosons) plays the 
role of a coolant [5|. Another interesting realization 
of boson-fermion mixture has been demonstrated in a 
two-component fermion system, where strongly bound 
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Cooper pairs correspond to bosons interacting with un- 
paired fermion atoms Starting from a microscopic 
description of a boson-fermion mixture we derive an ef- 
fective action for the order parameter (the condensate 
density) taking into account fermion-boson interaction. 
We show explicitly that the classical potential for the or- 
der parameter due to such interaction has two minima 
corresponding to the two phases of the system (mixed 
and phase separated) . The two minima are separated by 
the finite energy barrier, which points out that such a 
transition is indeed of the first order [7j . We then derive 
an expression for the nucleation (tunneling) rate of the 
critical droplet of the pure fermion phase near the phase 
transition line and near the line of absolute (spinodal) 
instability of the mixed phase. We also evaluate the role 
of dissipation [8] in the quantum nucleation process and 
find that near the line of the first order phase transition 
it changes the leading asymptotic behavior of the nucle- 
ation rate on the degree of supersaturation. 

We consider a Bose-Einstein condensate interacting 
with a single species of fermions (in the same spin 
state). Interactions in such a mixture are character- 
ized by two scattering lengths asB and ubf- Fermions 
and bosons interact through contact potential Asi?(5(r — 
r'), contributing term XbfP'4'p''Pf, where Xbf ~ 
2iTh^aBF{^/'mB -f l/m_F); y^e"^ and ij^F are bosonic and 
fermionic fields respectively. In addition, boson-boson 
interaction gives rise to another term \bbP^ I'i- in the 
Hamiltonian density, with \bb — ^T^h^aBB/fns- The 
direct coupling between fermions is negligible due to the 
exclusion principle in the s-scattering channel (p-wave 
scattering is usually small compared to s-wave fermion- 
boson and boson-boson interactions). For the purposes 
of present calculation we can neglect the spatial depen- 
dence of the trapping potential and assume that the local 
densities of fermions and bosons are set by the chemical 
potentials ^f and /is. Indeed, since the first order tran- 
sition occurs at finite coherence length (to be defined 
below), the shape of the trap potential should play lit- 
tle role in the dynamics of the phase transition as long 
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FIG. 1: (a) The phase diagram of the equihbrium 
boson-fermionic mixture. The bosonic density ub is in 
units of '^(a%g/a%p)(rn%/Tn%) (rnp/mB + and the 

fermionic density np is in units of '^{agg /a%p)(rn%/m%) 
X {mp/niB + The area where pure fermionic fraction 

can coexist with the mixture (F-M) is hatched. The mixture 
is unstable above the dashed hue. (b) Three possible curves 
for E{p) at equilibrium (sketch), (c) The sketch of E{p) for 
out-of-equilibrium mixture: above the lower phase separa- 
tion curve (dashed curve) and below the instability line (solid 
curve). The corresponding points are shown on the phase 
diagram (a). 



as the effective size of the trap is much greater than the 
coherence length. 

It is convenient to describe the system in terms of the 
bosonic field only, averaging e"^/*^^^ (H is the overall 
Hamiltonian) with respect to the fermionic field. Such 
averaging can be easily carried out within mean field, 
i.e., the Thomas-Fermi approximation so that 

the calculation reduces to the evaluation of the canonical 
partition function (or free energy) of the free fermions 
with effective chemical potential fxp — ^bfP- In the zero- 
temperature limit the effective classical potential density 
for the bosons is 
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where 9{x) is a step function. It is known that a fermion- 
boson mixture with pointlike interactions exhibits three 
different phases uniform mixture of bosons and 

fermions (M), pure fermion fraction coexisting with the 
mixture (F-M), and separated bosonic and fermionic frac- 
tions (F-B), see Fig[T^. Equation ([T]) is sufficient to ob- 
tain the entire structure of the phase diagram. Indeed, 
one can notice that E{p) has either one or two local min- 
ima separated by the barrier (see Fig. [T|d). The one at 
p — corresponds to pure fermionic phase. The other, 
p — Pa, characterizes the mixture if po < p-f/^bf, or 
pure bosons, if po > pf/^bf- At low densities only the 
mixture exists at equilibrium. For higher densities the 



mixture becomes metastable and, eventually, unstable 
above the absolute (spinodal) instability line. The phase 
transition line between M and F-M phases can be ob- 
tained from the condition E{0) = E{pq) in Eq. ^ in the 



parametric form: np = A ~ A ^[1 — ?/^]/4 



F — ^ y /"j 
where the densities are in the units of Fig. [1^, A = 

37r^/i'^ABB/(2mj?)^/^A^j^/i^''^, and parameter y is ob- 
tained from the equation 2 + 4:y + 6y^ + 3y^ — 5A{l + y)^ = 
0. Note that the solution with < y < 1 exists only for 
2/5 <A< 3/4. 

In what follows we consider the system with fixed 
global density of the bosons ub- Upon variation of 
fermion density in the metastable region, i.e., between 
points X and Y in Fig. [T^, the classical potential E{p) 
varies continuously between the two situations shown 
schematically in Fig.[T}:: with nearly equal minima in the 



vicinity of the phase transition point (at np 
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and vanishingly small barrier near the absolute instabil- 
ity at = Up = 8/27 (in dimensionless units). In the 
rest of this Letter we study kinetics of the system, i.e., 
the rate of formation of the stable pure fermion phase 
out of the metastable mixed phase in these two limiting 
cases. To do so we consider an effective Lagrangian den- 
sity of the bosons, which can be written in density-phase 
variables as 
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The first term in Eq. ^ is the Berry phase term, while 
the second and the third terms are the kinetic energy of 
the superfluid. The Thomas-Fermi approximation uti- 
lized in Eqs. (|ll2p implies that renormalization of the 
boson kinetic energy arising due to the non-locality of 
the fermionic response function is relatively small. A 
straightforward perturbative estimate (to the second or- 
der in Xbf) yields the gradient term, i.e. the correction 
to the Thomas-Fermi, - {m]ih%ppL]l'^ /h^kl,)(y pf . 
Comparing this term with the third term in Eq. Q we 
see that near the phase transition line it is smaller by 
factor ^ (kpl)^'^, where I is boson coherence length, 
I ~ o-BB/gB, and = u^bbPo is conventional boson 
gas parameter, which, in terms of the dimensionless ub 
(as in Fig.l) is 
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For gB ^ 0.1 and not too small np (e.g. for ub ^ 0.4), 
{kplY ^ 50, and thus the Thomas-Fermi approximation 
is well justified. 

The decay rate per unit volume (T/V) from a 
metastable state (at po) can be obtained by calculat- 
ing the classical action for the transition between states 
with p — po and p = in the imaginary time formal- 
ism by following the prescription of Ref. 0]. Namely, 
T /V ~ exp(— S*/?!,), where the action S = J dtdrL{p, (j>) 
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is evaluated over the classical (extremal) trajectory, de- 
fined by equations dS/6(j) = and 6S/Sp = 0. The first of 
these equations is a continuity equation, dtp+V {pu) = 
with u = hV(t>/mB- It can be easily solved in the spher- 
ically symmetric case yielding u = {r/r^p) dr'r' dtp. 
In terms of single variable p the action can be cast in the 
form 
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Evaluation of the extremum of S in the two cases of 
interest, i.e., near the line of the phase transition (at 
point X in Fig. [Iji) and near the absolute instability line 
(at point Y), has been formally carried out in Ref. 
In the former case the solution can be parameterized as 
p{r,t) « po^[^ ~ -R(^)] (the thin wall approximation), 
where R{t) is the the radius of the critical droplet and po 
is bosonic density of the mixed phase as before. Within 
such an approximation the action S can be formulated 
in terms of R{t) as 
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where the effective surface tension coeffi- 
cient is fj = ^/Wj2rriB J d\fp\/ E{p) and for 
IS.np/nPp <C 1 we have A/i — XBpKAnp, with 

K — \'^(n',?^\nn)'^/'^—nBn'p^^—n^pynBn'p^^, where 



ub and Up are dimensionless as in Fig. [1^ [ub varies 
between and 25/16). Evaluating the surface tension 
coefficient according to Eq. ([1]) and extremizing the 
action Sx, one obtains [11 1 



InFx/Fo 



11/2 



-0.0056- 



9B 



i^7/2 I Anp. 



7/2 



(5) 



Precise evaluation of coefhcient Fq lies out of the scope 
of present calculations. It can be estimated, however, 
as Tq/V ~ ojq/I^, where is an "attempt" frequency. 
From the uncertainty principle wo ^ h/2mBl^ and thus 
To/V - h/rriBl'^. 

As expected the tunneling exponent, i.e. the right- 
hand side of Eq. (O, is singular in the degree of metasta- 
bility Anp and diverges as Aup'^^^ Equation ^ 
also indicates that the rate of nucleation is exponentially 
small in the dilute limit, i.e., for qb <C 1. Note that ap- 
pearance of the dimensionless boson density in the 
numerator in Eq. ^ does not fix this problem, since the 
thin wall approximation (nucleation) requires sufficiently 
high energy barrier, e.g. Fig.[Tl^b,c), which is not the case 
when Ub ^ 1. However, due to the smallness of the nu- 
merical coefficient in the right-hand side of Eq. ([5]) one 
can hope that quantum nucleation is observable in suf- 
ficiently strongly coupled systems (which are presently 
realizable with the use of Feshbach resonance). Indeed, 



for gB ~ 0.1, nB ^ 0.4, and Anp/rip = 0.15, the coeffi- 
cient K 0.27 and the tunneling exponent is ~ —27. For 
the same parameters and Qbb 20 a.u., Tq/V can be 



estimated to be ~ 10^^ s ^pm ^ 



lO^^s ^pm which yields nucleation 



rate Tx/V ^ 1 s^^pm^^. 

Let us now evaluate the rate of transitions near the 
absolute instability of the mixed phase (at point Y in 
the phase diagram in Fig. [1^). Since the energy barrier 
in E{p) is now relatively small and E{po) E{0), see 
Fig- Hfc, one should pay attention only to the vicinity 
of the metastable minimum at po. Expanding action in 
Eq. ([3]) in Sp, p = po + 6p, and retaining the second and 
the third order terms (in Sp) in potential energy E{p) 
and only the second order terms in the kinetic energy 
(the third order kinetic terms contain gradients and thus 
are small compared to the the third order potential terms 
in the limit of diverging coherence length near the point 
of absolute instability) one obtains 
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The extremum of Sy is difficult to evaluate exactly. 
However, by introducing dimensionless variables x — 
ry/SniBPaa/h, r = Apoat/h, and p — 6pb/a, the action 
Sy can be rewritten as const x sy, where si'[p(x,r)] 
is a parameter-independent functional, the extremum of 
which is a c-number. Its value can be estimated by varia- 
tional ansatz p = —pQ exp (— ax^ — I3t'^), where a, (3 and 
Po are variational parameters. Upon a straightforward 
calculation one obtains 
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Again we see that tunneling exponent is controlled by 
the inverse boson gas parameter gB- The exponent van- 
ishes when fermion density np reaches value n^, where 
the effective energy barrier disappears, see Fig. [TJ^. In 
this respect Eq. ([7]) is similar to results on macroscopic 
quantum tunneling (MQT) in systems of trapped bosons 
with attractive interactions |12| . but with an important 
distinction: unlike in the latter case, the height of the 
potential barrier and thus the tunneling exponent for the 
critical droplet are controlled not by the total number of 
particles in the trap but the local densities. Therefore, 
for sufficiently large numbers of particles in the trap the 
tunneling exponent can be fine-tuned with a desired ac- 
curacy, which makes it possible to observe the MQT rate 
in a well controlled and predictable regime. 
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We now consider effects of dissipation on the kinetics 
of the first order phase transitions. The Thomas-Fermi 
approximation utihzed in the derivation of the effective 
potential E{p) impHes that fcrmions instantaneously ad- 
just to the local variation of boson density. Therefore it 
is instructive to proceed beyond this approximation and 
evaluate the effect of the excitation of particle-hole pairs 
in the fermionic subsystem on the dynamics of the tran- 
sition. This effect can be estimated by considering the 
second order frequency dependent correction (in boson 
density p) produced by the interaction XbfiPI^iPfP- The 
resulting contribution to the effective action can be cast 
in the form 



dq duj 
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where the oj and q independent term (omitted above) 
have already been included in Eq. ([T]), as can be verified 
by expanding E{p) to the order 0{p^). Again we con- 
sider the two cases: near the phase transition and near 
the absolute instability. Utilizing the thin wall approxi- 
mation, p(r,t) « po(^[f ~ ^6 obtain that near the 
line of the phase transition the action in Eq. ([5]) can be 
written in terms of the radius of the critical droplet as 
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The first two terms in the right-hand side of Eq. © arise 
due to the restructuring of the fermionic density of states 
inside the droplet in the course of its expansion, while the 
last term can be viewed as coupling between droplet 's 
surface and particle-hole excitation in Fermi sea. Exact 
evaluation of extremal action Q with the correction ^ 
is difficult and therefore we use the variational technique 
to evaluate the tunneling exponent. A natural ansatz is 
R{t) = i?oe~"* , where the turning point Rq and coeffi- 
cient a are variational parameters . We find that correc- 
tion due to ASx strongly alters the tunneling exponent 
in the region KAnF/nF<t:T]dis^ g^B^ni^ (mF/mB)^/^: 
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while in the opposite (nondissipative) limit Anp/rip 3> 
rjdis the tunneling exponent is given by Eq. We see 
from Eq. pO)) that the influence of dissipation on the 
dynamics of nucleation is significant for tub < mp- 

Near the absolute instability (for np^rip) the con- 
tribution of dissipation, e.g. Eq. ([H]), can be estimated 
by using the same variational ansatz as in derivation of 
Eq. ([7]) . We find that the tunneling exponent acquires an 



additional term ^ —{mpgBn^/mB)^^^ (the numerical 
coefficient in Eq. ([7]) does not change significantly when 
AS* is taken into account). This term is independent of 
fermion density and is again controlled by the mass ratio 
mp/niB- Therefore for not too high fermion/boson mass 
ratio, it is of the order g^^^ and thus dissipation does not 
significantly alter the dynamics of the phase transition 
(MQT) in the absolute instability region. 

Finally we estimate the crossover temperature between 
thermal and quantum regimes for the transition dynam- 
ics in the two cases of interest. The thermal regime be- 
comes effective when the energy gap between the exited 
and the ground state energies of the metastable mix- 
ture is of the order fc^T [l3|. This energy difference 
is HTq^ lT?/2mBl^- Near the phase transition curve it 
gives the relation g^^UB^ T/Tc, where Tc is the Bose- 
Einstein condensate transition temperature. Since we 
assume T<^ Tc and ub 1, the tunneling is clearly the 
dominant process of nucleation. Near the spinodal insta- 
bility we obtain g^^^ [l — np / n^p) jTc- In this case, as 
expected, thermal excitation energy bounds the coher- 
ence length l'^ {'^ — np/rip)''^/'^ for which the transition 
is dominated by MQT mechanism. 
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